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Recent work has found that the interplay between spin accumulation and Coulomb blockade in 
nanoparticle arrays results in peaky I-V and tunneling magnetoresistance (TMR) curves and in 
huge values of the TMR. We analyze how these effects are influenced by a polarization asymmetry 
of the electrodes, the dimensionality of the array, the temperature, resistance or charge disorder and 
long-range interactions. We show that the magnitude and voltage dependence of the TMR does 
not change with the dimensionality of the array or the presence of junction resistance disorder. A 
different polarization in the electrodes modifies the peak shape in the 1-V and TMR curves but not 
their order of magnitude. Increasing the temperature or length of the interaction reduces to some 
extent the size of the peaks, being the reduction due to long-range interactions smaller in longer 
arrays. Charge disorder should be avoided to observe large TMR values. 

PACS numbers: 



I. INTRODUCTION 

A large amount of effort has been devolted to under- 
standing and controlling the interplay between the elec- 
tronic current and magnet isirP^. Such an interplay is 
specially interesting in the presence of charging effect^. 
At low temperatures and voltages, the transport through 
a small metallic island is blocked due to the cost in energy 
to add an electron to the island. At zero temperature 
to allow current flow a minimum voltage, the threshold 
voltage Vth has to be applied to the electrodes. When 
the island is placed in between two ferromagnetic elec- 
trodes and the spin relaxation time is large, spin accu- 
mulation can appear at the island, i.e. a spin splitting of 
the chemical potential is createcP. Such spin accumula- 
tion is a non-equilibrium process which appears at finite 
voltages induced by the current flow. The spin accumu- 
lation happens to equalize the spin dependent tunneling 
rates to enter and exit the nanoparticle. 

When the polarization of source and drain electrodes is 
the same spin accumulation appears for antiparallel (AP) 
arrangement of the electrode magnetization, but not for 
parallel (P) orientation and induces a voltage dependent 
TMR^s, defined as 

TMR = ^^P-^P = Ip_jAP 
Rp Iap 

When the spin polarization of the electrodes is different 
spin accumulation also appears for P magnetizations and 
modifies the value of the magnetoresistance. With in- 
creasing temperature charging effects are reduced and the 
I-V curves and TMR evolve towards the non-interacting 
values. The specific device configuration determines the 
transport characteristics and even regions with weak neg- 
ative differential conductance or oscillations in the TMR 
can be founcPtiSI. 

Chargin g effe cts are known to be enhanced in nanopar- 
ticle arrayd^'i^lti^. The threshold voltage increases, being 
proportional to the number of islands when interactions 
are restricted to charge s in t he same conductor or when 
there is charge disordei^^E^. For onsite interactions all 



the voltage drop happens at the contact junctions, be- 
tween islands and electrodes. The large threshold volt- 
age is a consequence of the charge gradient which has to 
be created at the inner junctions, between the nanopar- 
ticles, to allow current. Accumulation of charges is the 
only way to create a potential drop at these junctions 
in clean arrays. Increasing the bias voltage does not in- 
fluence this potential drop if the charge state does not 
change. This changes for long-range interactions. In this 
case there are junction dependent voltage dropJ^. On 
the other hand charge disorder induces voltage indepen- 
dent potential steps at the tunnel junctions and modifies 
the threshold voltage. 

Recently we showed that when a nanoparticle array 
is placed between two ferromagnetic electrodes, the in- 
terplay between charging effects and spin accumulation 
has a dramatic effect in the current and TMR^^. Non- 
homogeneous spin accumulation appears for P orienta- 
tion of the electrodes. The inhomogeneity of the spin po- 
tential favors a reduction of the threshold voltage thanks 
to the spin-dependent potential drops created at the in- 
ner junctions. Because of the relation between spin ac- 
cumulation and current, large oscillations appear in the 
I-V curves. The spin accumulation decreases when a 
new conduction channel is open^^ producing a reduction 
in the current with increasing voltage, and oscillations 
appeai^^. The anomalous behavior is found below the 
non-magnetic threshold voltage. Spin accumulation is 
also present for AP arrangement of the magnetizations, 
but it is fairly homogeneous, so the threshold voltage 
and the overall I-V curves do not change much when 
compared with the non-magnetic results. The different 
threshold voltage for P and AP configurations results 
in huge values of the TMR, which are strongly voltage- 
dependent. A different behavior in P and AP configura- 
tions is also observed in the relative polarization of the 
current (/^ — ■ When the polarization of source and 
drain electrodes is equal, the current is not spin polar- 
ized for AP arrangement as it also happens in the single 
particle case. On the contrary for P configuration the 
spin polarization of the current presents clear oscillations 
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which correlate with the ones found in the currentP^. 

The fact that these huge values of the TMR are associ- 
ated with the (in absence of spin-accumulation) vanish- 
ing differences in energies for tunneling through the inner 
junctions suggests that this effect could be extremely sen- 
sitive to modifications of the array characteristics, such 
as length of the interaction, presence of disorder, asym- 
metry in the electrodes polarization, among others or to 
external parameters such as temperature. In this pa- 
per we analyze such sensitivity. We have found that the 
magnitude and voltage dependence of the TMR do not 
change with the dimensionality of the array or the pres- 
ence of junction resistance disorder. Polarization asym- 
metry modifies the peak shape in the I-V and TMR 
curves but not their order of magnitude. The size of 
the peaks is reduced to some extent with increasing the 
temperature or length of the interaction but the effect 
of long-range interactions is reduced in long arrays. On 
the other hand charge disorder is harmful and should be 
avoided in order to observe these effects. 



II. THE SYSTEM 

We consider an array of metallic nanoparticles placed 
in between two ferromagnetic electrodes. Here the word 
nanoparticle, also called island, refers to confined metallic 
regions separated from each other and from the electrodes 
by tunnel junctions, not necessarily sphere like metallic 
particles. They can be nanostructures created by other 
methods such as lithography or epitaxy, for example. 

Due to confinement to add an electron to the system 
costs a charging energy Ef.- We restrict ourselves to the 
classical Coulomb blockade regime with 5 <C KbT <C 
Ec with 5 the level spacing, T the temperature and Kb 
the Boltzmann constant. The arrays are one-dimensional 
(ID) and contain N nanoparticles, see sketch at the top 
of Fig. [ij except in section IV where the effect of the 
dimensionality is studied. In this section two-dimensional 
(2D) arrays with MxN particles have AI rows containing 
each of them N nanoparticles, see sketch in Fig. [3j 

The tunnel resistances which separate the nanoparti- 
cles from other nanoparticles are — 2Rt for both spin 
orientations through all the paper, except in section |VI| 
where they are allowed to vary randomly between two 
values to analyze the effect of resistance disorder in the 
TMR. In equilibrium the nanoparticles are non magnetic 
but the source and drain electrodes are ferromagnetic 
with spin polarization pi and p2. The polarization of 
the electrodes enters via the tunneling resistance which 
depends on spin. At the contact junctions between the 
array and the electrodes = 2Rt{1 ±pi)~^ with a the 
spin index. Plus (minor) signs are assigned to major- 
ity (minority) spin carriers and i — 1,2 to source and 
drain electrodes. Here 2Rt is to be substituted by the 
assigned resistance value if disorder is included. We con- 
sider Pi = P2 = P, except in section |III| where the effect 
of the polarization asymmetry of the electrodes is stud- 
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FIG. 1: Top: sketch of the nanoparticle array under consider- 
ation. The arrow at the drain refers to the P (light) and AP 
(dark) configurations. (a)TMR as a function of the voltages 
corresponding to an A^ = 10 clean array with source polar- 
ization pi =0.7 and different drain polarizations p2. (b) and 
(c) I-V curves for the same arrays in (a) for antiparallel and 
parallel orientation respectively. 



ied. Transport is treated at the sequential tunneling level 
with tunneling rate^^fil 



r.(Ai;„) 



^E„ 



R, expi^E./KeT) - 1 



(2) 



where the electronic charge has been taken equal to unity. 
A£^cr is the change in energy of electrons with spin a due 
to the tunneling process and can be written as 



{4>P - 4>a) 



(3) 



where indices a and /? label the two conductors involved 
in the tunneling. 0^ and 4>f3 refer to the potential at the 
conductor from which the electron leaves and at which 
the electron arrives, respectively. E^~i^ is the cost in 
energy to transfer an electron between a and f3 assuming 
that the array is clean, at zero potential, uncharged and 
in equilibrium and equals 



E°' + E^ 



(4) 



For the islands E" — Ec, the unit of energy, while at 
the electrodes E^ <^ Ec and is neglected in the discus- 
sion, is the matrix which describes the interaction 
between electrons. The diagonal elements = 2E". 
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In Eq. Q C~p measures the attraction of the electron 
and hole created in the conductors due to the tunneling 
process and is finite only when the range of the inter- 
action does not vanish, long-range limit. Below we con- 
centrate in the short-range interaction limit = Sap 
in which interactions between electrons are restricted to 
those charges placed in the same conductor. In sec |VIII| 
we use = C~^el"~'^l/°" to analyze how long-range 
interactions influence the TMR. An exponential decay is 
a good approximation when the islands are capacitively 
coupled to the neighboring ones or for screened interac- 
tions. Similar approximations h ave bee n frequently used 
in studies of nanoparticle arrayJ^^l^^. Here ao can be 
understood as an effective interaction or screening length 
and is measured in units of the nanoparticles center-to- 
center distance. 

At the source and drain electrode the potential are 
(!> source = V /2 and (t>drain = -V/2 respectively. The 
potential at the islands can be decomposed in different 
contributions 
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Here 



Na a-S with Na a the number of excess elec- 



trons with spin a at island a. Spin accumulation at the 
islands is defined as — ^"'^ = [Na^a — Na-a)S and 
it can be finite under a non-equilibrium current. The 
charge contribution to the potential depends on the total 
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essentially plus a small modification due to 
the proximity of the electrodes at a fixed potentiaP^. 
In the short-range limit 0^'' = 2EcQa- 0q°' is the po- 
larization potential at the island, induced by the elec- 
trodes at finite bias. In the short-range case it van- 
ishes. This means that there is no polarization poten- 
tial drop at the junctions situated between two islands; 
the polarization potential drop is only finite at the con- 
tact junctions between an island and an electrode. In 
the presence of long-range interactions 0^°' is propor- 
tional to the bias voltage. For exponentially decaying 
interactions 0^°' = ^{6'°'^''° - e-l^+i-"l/°o). The po- 
larization potential drop at the inner junction increases 
with increasing the length of the interactions^^. The last 
term in Eq. ^ c/)'^^ is a contribution associated to the 
presence of random charges in the substrate or array sur- 
roundings. It is finite at the islands in charge disordered 
arrays (sec. VII) and zero in clean arrays, even when 
there is resistance disorder. If interactions between the 
charges are short range (C~^ = 5a, p), the set of disorder 
potentials {(j)'^'^}, once the screening of the potential due 
to the mobile charges is taken into account, is uniformly 
distributed in the interval —Ec < 0*'' < Ec- In the 
presence of long-r ange interactions the screened disorder 
is correlate d^^l^^l Here we only consider the presence of 
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FIG. 2: (a) and (b) Spin accumulation at the islands as a 
function of the position for the array parameters in Fig. [l] 
and V = 16i?"' respectively for parallel and antiparallel ar- 
rangement of the electrode magnetizations, (c) and (d) Spin 
polarization of the current as a function of the voltage for the 
same arrays in (a) and (b). 



charge disorder 0'^*'' in the case of short-range interac- 
tions and the first situation applies. 

Except otherwise indicated we consider a ID array 
with short-range interactions, homogeneous resistances, 
absence of charge disorder and take p = 0.7, ksT — 10~^ 
and 5 = 10~^ with the energies given in units of Ec- The 
value of 5 does not influence the results as soon as 5 ^ T. 

The current is calculated by means of a Monte Carlo 
simulation described elsewher^^^f^^HH Simulations be- 
come increasingly costly with increasing temperature. 
When discussing the effect of temperature we restrict 
ourselves to small arrays and temperatures smaller than 

QMEr. 



III. ELECTRODE POLARIZATION 
ASYMMETRY 

In the case of a single island the polarization asymme- 
try ipi 7^ P2) induces spin accumulation for P orienta- 
tion, which is absent if pi — p2. While the dependence of 
the TMR on voltage does not change much qualitatively, 
its magnitude depends on the values of pi and p2 , satu- 
rating at high voltages to 4pip2/(4 — {pi +^2)^)- In ar- 
rays, to have different spin polarization in the electrodes 
breaks inversion symmetry and is expected to modify the 
spin potential profile created along the array. In particu- 
lar, it can create a spin potential gradient not only when 
the electrodes polarizations are parallel, but also when 
they are antiparallel. The presence of such a gradient 
could modify notably the I-V curves of the AP configu- 
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ration, allowing current flow below the metallic thresh- 
old and reducing the magnetoresistance in the voltage 
regime, where the peaks are observed. 

Fig. [ija) shows that the peaks in the TMR below the 
non-magnetic threshold {Vr — l8Ec for N = 10) do 
not disappear when spin asymmetry is present. How- 
ever, their shape is strongly modified. The values of the 
TMR have the same order of magnitude as observed when 
Pi = P2- Whether they are larger or smaller depends in a 
complicated way on the polarizations and on the voltage. 
Below Vt and for fixed pi the current found for parallel 
arrangement increases with increasing p2 without alter- 
ing much the shape of the I-V curve, see Fig.[l](c). With 
AP orientation peaks are observed in the I-V only when 
Pi 7^ P2 J see Fig. [l|b) . Their shape is more irregular and 
their height notably reduced as compared to those found 
for P-orientation. The change in the shape of the TMR 
curve is mostly due to the change in the current for AP 
arrangement. 

The differences in the current when compared with the 
Pi — P2 case originate in the different spin gradient cre- 
ated through the array by the spin accumulation. The 
change in spin accumulation between the first and last 
island of the array is much larger when the electrode 
magnetizations are P than when they are AP, as shown 
in Fig. [2] (a) and (b). For P-orientation when pi — p2 
the sign of the spin accumulation changes at the center 
of the array. With spin asymmetry this change of sign 
can happen or not and if it does, the change of sign will 
not happen at the array center in a generic case. For 
P arrangement and spin asymmetry the sign of the spin 
gradient does not depend on the asymmetry, but it is 
controlled by the orientation of the electrode magneti- 
zation. With AP orientation the spin accumulation can 
increase or decrease along the array depending on p2 be- 
ing larger or smaller than pi . At least for the cases that 
we have analyzed we have never found a change of sign 
in the spin accumulation. 

The polarization asymmetry also modifies how much 
spin polarized it is the current. As it also happens in the 
single island case, the current in the AP-arrangement can 
be spin polarized only if pi ^ P2 , with the polarization 
sign depending on p2 being smaller or larger than pi , see 
Fig. [2|^d) . For P-orientation it depends on the values of 
the electrode polarizations but does not change sign, as 
seen in Fig. 2(c). At high voltages the spin polarization 
saturates to 
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FIG. 3: (a)I-V curves in one and two dimensional arrays with 
parallel orientation: ID (N=5), 2D(N=3x5) and p=0.7. The 
I-V curve for the 2D system has been divided by 3, the number 
of rows, (b) TMR as a function of the voltage for the arrays 
in (a). 



trodes from N ^ 1.2toV = 3or larger has an enormous 
impact on the I-V curves and on the TMR. Thus it is sen- 
sible to ask what does it happens if we increase the width 
of the array to become two-dimensional. For simplicity 
we assume an square-lattice array, see inset in Fig. [3] 
The current through metallic two-dimensional clean ar- 
rays (no charge disorder) is reasonably well described by 
assuming that the different rows conduct in parallel^. If 
this is true also in the presence of spin accumulation no 
much change would be expected in the magnetoresistance 
because the current for both P and AP would be given by 
that of a one dimensional array multiplied by the number 
of rows. The validity of this approximation is confirmed 
in Fig. 3 where the TMR and the I-V curve for paral- 
lel arrangement for a ID N=5 array and for a 2D 3x5 
one (three rows of 5-particles length) are compared. The 
agreement between the I-V curves is very good, once the 
2D-array I-V has been divided by the number of rows. 
Minor differences are found around the metallic thresh- 
old {Vt = 10 for N = 5). As expected the TMR of the 
2D arrays follows very well the one of the ID array. A 
small difference is only seen at the top of the peak. The 
TMR is more sensitive to small changes in the current 
of the AP-configuration (not shown), as this one is much 
smaller and enters into the denominator. 



V. TEMPERATURE 



IV. DIMENSIONALITY OF THE ARRAY 

As shown irfiSl to increase the length of a one- 
dimensional array placed between ferromagnetic elec- 



Thc high values of the TMR at the peaks are expected 
to be quite sensitive to temperature. These values are as- 
sociated with the small current which flows through the 
array below the non-magnetic threshold when the mag- 
netizations are AP. Such a current extrapolates to zero at 
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FIG. 4: (a) TMR as a function of temperature for a clean 
iV = 5 array and a voltage V = 8.4i5c. Top inset: TMR as 
a function of temperature for a single island and V = 2Ec. 
It shows a much weaker temperature dependence, what is 
emphasized in the lower inset where the dependence of the 
relative values of the TMR as a function of temperature are 
represented, (b) and (c) I-V curves for an A'' = 5 array with 
AP and P orientations, respectively, for two different temper- 
atures. The effect of temperature is much stronger for the AP 
orientation, as expected. 



zero temperatures. Such vanishing current at low tem- 
peratures is associated with tunnchng processes at the 
inner junctions with vanishing cost in energy. Increasing 
the temperature allows the tunneling event to happen. 
From Eq. [2] the tunneling rate of processes with zero en- 
ergy cost increases linearly with the temperature. On a 
first approximation the current for AP orientation is ex- 
pected to increases linearly with the temperature, while 
the current through an array with P magnetizations will 
be much less sensitive because the tunneling processes at 
the inner junctions have a finite energy gain provided by 
the spin potential gradient. Such dependencies would re- 
sult in a TMR inversely proportional to the temperature. 
Confirmation of this dependence and on the different ef- 
fect of the temperature for P and AP orientations is seen 
in Fig.|4l 

Note that the strong dependence of the TMR on the 
temperature is not a simple consequence of the weaken- 
ing of Coulomb blockade. As shown in the top inset of 
Fig. |4] (a) the temperature dependence of the TMR corre- 
sponding to the single island case = 1 is much weaker 



at KbT <^ Ec, because the tunneling processes involved 
are thermally activated, contrary to the zero energy cost 
relevant at the inner junctions of a long array with AP 
orientation. The different temperature dependence of the 
TMR for a single island and a long array is emphasized 
in the lower inset of the same figure. 



VI. RESISTANCE DISORDER 

The strong exponential dependence of the tunneling 
probability on the tunneling junction width makes that 
small differences on the distance between the nanopar- 
ticles result in large differences between the junction re- 
sistances. As a consequence resistance disorder is ex- 
pected to play an important role in many devices based 
on nanoparticle arrays. In this section we show that the 
appearance of peaks in the I-V curves and the large val- 
ues of the TMR are robust against resistance disorder. 

Fig. [5] (a) shows the I-V curves for P and AP config- 
urantions corresponding to an TV = 20 array with resis- 
tance disorder, being the junction resistances randomly 
assigned and varying between (10 — 21)Rt- Due to the 
different effect of the change in resistances for different 
voltages, these I-V curves are not compared with those 
found in clean arrays. Qualitatively, the main features 
displayed are the same in both cases. The I-V curve cor- 
responding to AP arrangement shows a threshold volt- 
age equal to the one found with non-magnetic electrodes, 
while the one for P-orientation is reduced. Peaks in the 
current are found at low voltages, below Vr, the non- 
magnetic threshold, when the electrode magnetizations 
are parallel. 

The TMR corresponding to these I-V curves is shown 
in Fig. [5] (b) where it is compared with the one found 
in the absence of resistance disorder. In the peaky re- 
gion of voltages no differences are found between the 
TMR of both systems. The high values of the magne- 
toresistance are extremely robust against resistance dis- 
order. As could be expected from the shape of the I-V 
curves m Fig. [5] (a) the spin accumulation, shown in Fig.js] 

(c) and (d) is quite homogeneous through an array with 
AP orientation and decreases from left to right with P- 
orientation. As it happened when all the resistances are 
equal it changes sign at the center of the array, see Fig. [3] 

(d) . For the voltages shown the spin accumulation dif- 
fers very little when compared to the disorder free case. 
Interestingly, and opposite to what happened in uniform 
arrays, in the AP configuration, at high voltages an spin 
gradient is formed through the array what produces slight 
differences between the TMR of clean and disordered ar- 
rays at these voltages, not shown. 

The spin polarization of the current is shown in Figs. [5] 

(e) and (f) for both magnetic orientations. In the P- 
arrangement it equals the one found when there is no 
resistance disorder. The situation changes with AP ori- 
entation. At high voltages a finite spin polarization is 
found in the disordered case, while it vanishes in the 
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FIG. 5: (a) I-V curves for a A'^ = 20 array with the junction 
resistances randomly assigned varying between (10 — 22) Rt, 
for both magnetic configurations and p — 0.7. (b) TMR as a 
function of the voltage for the same array that in (a) compared 
with the one corresponding to a resistance disorder free array, 
(c) and (d) Spin accumulation at the islands for V = 38.5Ec 
with antiparallel orientation and V — 20Ec and parallel orien- 
tation, respectively for disorder free and resistance disordered 
arrays, (e) and (f) Spin polarization of the current for antipar- 
allel (d) and parallel (e) arrangements for disorder free and 
resistance disordered arrays. The curves in all the graphs cor- 
respond to the same N = 20 arrays with or without resistance 
disorder. 



VII. CHARGE DISORDER 

Besides resistance disorder, charge disorder can be im- 
portant in the devices. This type of disorder is un- 
avoidable in many systems and specially important in 
self- assembled nanoparticle arrays. It arises as a conse- 
quence of nearby charge impurities which create random 
potentials in the nanoparticles. These random potentials 
can be described in terms of random background charges 
spread in the interval (—1/2, 1/2) in each nanoparticle. 

Charge disorder has a strong impact on the I-V curves 
of nanoparticle arrays even when the electrodes are non- 
magnetic. It modifies the threshold voltage, which now 
depends on the particular disorder configuration. In 
the limit of short range interactions discussed here, the 
threshold voltage, on average, is reduced to the half of the 
value found in a clean-array. Only half of the junctions 
(upward steps in potential) prevent the flow of charge. 
The voltages at which steps are found in the Coulomb 
staircase regime are also disorder dependent. 

Fig. |6] (a) compares the I-V curves of clean and disor- 
dered arrays when the electrodes are magnetic and their 
polarizations are parallel or antiparallel. The I-V curves 
of the disordered arrays strongly differ when compared 
with the clean case. Most significatively, the strong de- 
pendence of the current on the magnetic orientation of 
the electrodes has disappeared. This is clearly seen in 
Fig. [6] (b) which shows a very small magnetoresistance. 
This figure is to be compared with Fig. [5](b). Random 
charge disorder has reduced the TMR up to four orders 
of magnitude. 

The large values of TMR found in clean arrays are a 
consequence of the spin potential gradient generated in 
the array present (absent) when the electrodes magneti- 
zations are parallel (antiparallel). This potential gradi- 
ent allows current to flow through the inner junctions. 
In the absence of such gradient, at zero temperature, the 
current vanishes at low voltages because, while tunnel- 
ing does not cost energy it also does not gain it. The 
suppression of the TMR is not due to a reduction of the 
spin gradient. With charge disorder spin gradients are 
present, being even larger than in the clean case, see 
insets in Fig. [6](b). They are finite, but much smaller 
with antiparallel configuration. However, contrary to the 
clean case, with disorder, tunneling through the junctions 
which prevent current (upward steps) costs energy. The 
change in potential due to spin accumulation is too small 
to overcome such energy cost or to modify significatively 
the tunneling rates through the junctions with downward 
steps in disorder potential. 



disorder-free array. The sign and dependence on voltage 
change for different realizations of the resistance disor- 
der. The finite spin polarization found at low voltages 
is consequence of the finite temperature which allows a 
small current fiow below the zero temperature threshold. 



VIII. LONG-RANGE INTERACTIONS 

So far we have studied onsite interactions in which the 
interactions between the electrons are limited to those 
charges placed in the same conductor. This case is spe- 
cial because all the voltage drop happens at the contact 
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FIG. 6: (a) I-V curves for A'^ = 20 clean and charge disordered 
arrays, with source and drain polarization p=0.7 with parallel 
and antiparallel configurations. (b)Main figure: TMR as a 
function of the voltage for the charge disordered array used 
in (a). Insets. Comparison of the spin accumulation gradient 
created in the clean and charge disordered arrays in (a) at 1/ = 
20-Ec with antiparallel (left) and parallel (right) orientations. 



junction. There is no polarization potential at the is- 
lands. This makes that at zero temperature and below 
the threshold corresponding to non-magnetic electrodes 
the current is maintained only by the spin accumulation. 
The potential drop at the inner junctions increases with 
increasing range of the interactions, as can be observed in 
Fig. [t] (a), where an exponentially decaying interaction, 
typical of capacitively coupled islands has been used, oq 
can be understood as a screening length. 

For finite range of the interactions and in the absence 
of charge disorder, once a charge has entered the array 
it is able to flow, even in the absence of charge gradient. 
The threshold voltage is strongly reduced and given by 
the voltage which a llows an electron to enter the array 
or a hole to leave it! ^'' l ^^ l The role played by the spin 
accumulation on the transport is expected to decrease if 
the interactions are long range. The I-V curves corre- 
sponding to P and AP orientation, are shown in Fig. [7] 
(b) and (c) for the values of ao used in (a) . The thresh- 
old voltage is similar for P and AP configurations, even 
if at voltages close to threshold the current is still highly 
suppressed, specially for the smaller values of oq. The 
potential drop at the junctions, and correspondingly the 
I-V curves shape are influenced both by the long-range of 
the interactions which affects the step width and height 
and by the spin accumulation responsible of the oscilla- 
tions. For small ao oscillations in the current are observed 
when the electrode magnetizations are parallel, but they 
disappear for large oq. 

The effect of increasing the range of the interaction on 
the TMR is plotted in Fig.jSj As expected the amplitude 
of the oscillations in the TMR is clearly reduced as oq in- 
creases. Note the different scale of the TMR axis in the 
plots. For TV = 10 and = 1 the TMR almost vanishes. 
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FIG. 7: (a) Polarization potential at the islands and elec- 
trodes for an N=10 array sXV = 16-Ec and several values of 
ao, the range of the interactions, ao = refers to the onsite 
interaction limit discussed in previous sections, (b) and (c) 
I-V curves corresponding the same values of ao and N in (a) 
for parallel and antiparallel orientation, respectively. 



On the other hand, as seen in Fig. [s] (c) the reduction of 
the TMR for a given interaction length depends strongly 
on the array length. As longer it is the array the reduc- 
tion in the magnitude of the oscillations is weaker. This 
dependence appears because given an interaction length 
the polarization potential drops at the junctions between 
particles is larger for smaller iV. 



IX. SUMMARY AND DISCUSSION 

In summary, we have seen that the oscillations in the 
TMR and in the I-V curves are robust to most of the 
modifications that can appear in an experiment, when 
compared to the ideal case studied irfiSl. In particular 
if source and drain have different spin polarizations the 
magnitude of the oscillations does not change much. The 
main effect induced by the polarization symmetry is a 
change in the peak shape and the appearance of peaks 
when the magnetizations are antiparallel. To go from a 
ID to a 2D-square lattice array does not affect the TMR 
to a first approximation, while the current is increased 
following the number of rows in the 2D array. A similar 
situation is found in arrays with resistance disorder. 

The amplitude of the oscillations and the TMR de- 
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FIG. 8: (a) to (d) TMR as a function of voltage for an AT = 10 
array for ao = 0,0.3,0.5 and 1, respectively. In (a) and (c) 
the TMR of an A'' = 5 array with the same ao is included. 
The reduction of the oscillations amplitude with ao is faster 
in shorter arrays. 



creases with temperature so, the temperature should be 
kept as small as possible in an experiment as compared 



to the charging energy. Both large and small values 
of KbT are convenient. 

In many arrays the long-range part of the interaction 
is screened by mobile charges in neigbour conductors. If 
this does not happen the TMR can be suppressed with 
respect to the onsite interactions case. This effect is less 
important in long arrays. Thus, if the length of the inter- 
action has to be taken into account the oscillations and 
large values of the TMR will be better observed in long 
arrays. 

The most harmful effect is charge disorder. In some 
systems like self-assembled arrays this kind of disorder 
seems unavoidable because of charges quenched on the 
substracted. However, in other devices like epitaxially 
grown pillars the disorder will be probably less impor- 
tant, what makes them a better candidate to observe 
these effects. 
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